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In this report we discuss the organization of different levels of nature and the cor- 



\q . Abstract 

in 
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ly-) , responding space-time structures by the consideration of a particular problem of time 

' irreversibility. The fundamental time irreversibility problem consists in the following: 

how to reconcile the time-reversible microscopic dynamics and the irreversible macro- 
scopic one. 

The recently proposed functional formulation of mechanics is aimed to solve this 
problem. The basic concept of this formulation is not a material point and a trajectory, 
^ , like in the traditional formulation of mechanics, but a probability density function. 

Even if we deal with a single particle (not with an ensemble of particles), we describe 
its state as a probability density function. We justify this approach using measurement 
theory. 

A particular problem in the framework of the irreversibility problem is the deriva- 
tion of the Boltzmann kinetic equation from the equations of microscopic dynamics. We 
propose a procedure for obtaining the Boltzmann equation from the Liouville equa- 
tion based on the BBGKY hierarchy, the recently proposed functional formulation 
of classical mechanics, and the distinguishing between two scales of space-time, i.e., 
macro- and microscale. The notion of a space-time structure is introduced. It takes 
into account not only the space-time itself (i.e., a pseudo-Riemannian manifold), but 
also a characteristic length and time. The space-time structures form a hierarchy in 
sense that the initial values for the processes on the microscopic space-time structure 
(interactions of the particles) are assigned from the processes on the macroscopic one 
(kinetic phenomena). 
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1 Introduction. Irreversibility problem 



The laws of physics are formulated in terms of differential equations. Every physical law is a 
mathematical model of the corresponding level of nature. For the last centuries theoretical 
and mathematical physics made the great success in the revelation of these laws as well as in 
the consequent technical progress. However, in order to have a holistic picture of the world, 
we should not only describe its separate levels, but also understand, how these descriptions 
relate to each other. 

One conception is reductionism. It claims that the properties of a system can, in principle, 
be derived from the properties of its parts. For example, the dynamics of a gas can be 
derived from the dynamics of its atoms or molecules. But it seems that the reductionism 
does not work: a system as a whole has some properties which cannot be understood from 
the studying of its parts. For example, the temperature characterise the gas a whole and 
cannot be understood from the level of separate atoms or molecules. 

Another property of this type is the property of irreversibility of macroscopic systems in 
time. Laws of the microscopic motion, given by the Newton, Hamilton, Schrodinger equa- 
tions, quantum field theory equations etc., are time-symmetric, i.e., all physical processes 
can go to both directions of time. On the contrary, the laws of macroscopic motion, given by 
the Boltzmann, Navier-Stocks, other kinetic and hydrodynamic equations, the second law of 
thermodynamics etc., have a distinguished direction of time. 

Let us consider this problem more mathematically. The Newtonian formulation of clas- 
sical mechanics based on the notions of point and trajectory in the phase space. The fun- 
damental law of dynamics in this formulation is given by the Hamilton system of equations. 
Consider the system of N identical particles. Then the Hamilton equations have the following 
form: 



Here = qi{t) gGc M 3 is the position of the z-th particle at time t, pi = Pi(t) G M 3 is the 
momentum of the z-th particle at time t, H = H(x\, . . . , xn) is a Hamiltonian, Xi = (qi,Pi). 
The following proposition holds: 

Proposition 1. Let (q(t),p(t)) be a solution of the Hamiltonian system of equations 



OP, q(t) = ( qi (t),...,q N (t)), p(t) = ( Pl (t),...,p N (t)). Let H(q,p) = H{q,-p). Then 



(q(t),p(t)) = (q(—t),—p(—t)) is also a solution of the Hamiltonian system. 

The proof is straightforward. This proposition is a mathematical expression of the time- 
reversibility of the classical dynamics of particles, since the most of the fundamental Hamil- 
tonians satisfy the condition H(q,p) = H(q, —p). 

We suppose that G is a bounded region in M 3 with the volume V and a smooth boundary 
dG. We will consider the Hamiltonians of the form 



OH 



dH 



1,2,. ..,N. 







i>j 
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Here $(r) is an interaction potential of the particles, U(q) is an external potential, m > 
is the mass of a particle. We assume that $(r) is continuously differentiable, monotonically 
decrease (which corresponds to the repelling force) on (0,+oo), $(r) — > as r — > +00, 
$(r) — > +00 as r — > 0. £/(g) is continuously differentiable in G, U(q) — > +00 as g — > dG 
(i.e., the external potential does not allow the particles to leave the domain G). fi > is a 
dimensionless scale parameter which we will use later. 

Consider now a gas from the point of view of kinetic theory. In this case the gas is 
described in terms of a single-particle density function f(q,p, t). Let it be normalized on the 
volume V (i.e., f(q,p,t)/V is a probability density function), q e G, p G M 3 . One of the 
fundamental equations of the kinetic theory is the Boltzmann equation: 

?l = _L?L + E.?l + stf, (3) 

dt m dq dq dp 
f \P — Pi\ 

St/ = n / [f{q,P,t)f{q,Pi,t)- f(q,p,t)f(q,p 1 ,t)]dadp 1 , 

</R 2 xR 3 771 

where n > is the concentration of the particles (the number of particles in the unit of 
volume), da = rdrd(p. Also p' and p\ are the momenta that will have two particles long 
after the collision, if they had the momenta p and p\ long before the collision with the 
impact parameter of the collision (scattering) r and the polar angle <p (so, (r, <p) e M 2 are 
polar coordinates on the plane perpendicular to the relative velocity vector (p — pi)/m). 
Thus, p' = p'(p,pi,r,ip), p[ = p[(p,pi,r,ip), the dependence is defined by the two-particle 
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Hamiltonian without the external potential H2 = + + $( ~ 92 )• The expression St / 
is called the collision integral. 

This is an important nonlinear equation which describes the relaxation of the function / 
to the Maxwell distribution. One of the properties which can be easily proved is the so-called 
Boltzmann if-theorem: 

Proposition 2 (if-theorem). Let f(q,p,t) be a solution of the Boltzmann equation (Ej) and 
the quantity 

f(q,P,t)kif(q,p,t)dqdp, 

n<mt&n'rnl r*r\nn n mime c ) T'h&m !. 

dt 



Qy = G x M 3 , be well-defined (i.e., the integral converges). Then ^jf- > 0. 



This proposition states the entropy growth (the quantity S(t) = —H(t) can be regarded 
as the entropy of the gas) and, hence, the irreversible character of the gas dynamics. 

Thus, we have obtained two contradictory conclusions: if we consider the g clS clS 8b whole 
(in terms of the single-particle distribution and the Boltzmann equation), the dynamics is 
irreversible in time, while if we consider the gas as a system of a finite number N of particles 
(in terms of the point in the iV-particles' phase space Qy and the Hamiltonian equations), 
the dynamics is reversible. 

It seems that the reductionism does not work, the time irreversibility is a property of 
macrosystems which cannot be reduced to the microscopic level. So, the problem about 
another type of relation between the different levels of description arises. 
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2 Functional formulation of classical mechanics and 



measurement theory 

Recently, a new formulation of classical mechanics called functional formulation was proposed 
by I. V. Volovich [TJ [2] and is aimed to solve the irreversibility problem. This formulation 
is based not on the notions of a material point and a trajectory, but on the notion of a 
probability density function for a finite number of particles. Consider again a system of TV 
identical particles. In the functional formulation the state of this system at time t is given 
by a probability density function p(x\, . . . ,XN,t). The fundamental law of motion in this 
formulation is given by the Liouville equation 



where {•, •} are the Poisson brackets. 

Formally, the Liouville equation is time-symmetric. But the solutions of the Liouville 
equation have the property of derealization. For example, if we have one particle in a 
box with hard walls (i.e., N — 1, U(q) — if q G G, U(q) — +oo if q G), then the 
spatial probability distribution u(q,t) = J R p(q,p,t)dp tends to the uniform distribution in 
G. So, initially spatially well localized particle diffuses and fill the whole of the aivailablle 
volume. This derealization property can be regarded as the irreversibility. So, the point of 
the functional formulation of mechanics is introducing the irreversibility on the microscopic 
level by the assignment of the holistic statistical properties to a single particle. 

The argument of considering the probability density function rather than the material 
point as a fundamental object of mechanics is given by the measurement theory. Since the 
real (more precisely, irrational) numbers are unobservable, the notions of a material point 
and a trajectory, based on the concept of a real number, have not a direct physical meaning. 
The result of every measurement is a rational number (it is worthwhile to note that the 
justification of the p-adic and adelic mathematical physics also starts from this fact [3]). 
Moreover, every measurement has an error. So, we never know a state of the system as a 
point in the phase space. Every experimenter knows a state as a rational number and an 
error represented as a confidence interval and a confidence probability. This situation is 
naturally described by a probability density function. 

For example, let an experimenter performs some measurements of some observable X and 
represent the results of the measurements in the form 



Here x is an "actual" value of X, i.e., he do not know the value of X, but we know that it lies 
in the confidence interval (x — A 0) 997, x + A 0j997 ) with the probability 0,997. Alternatively, 
he can represent the result of the same measurements in the same form with some another 
confidence probability 7 and the corresponding A 7 (of course, A 7 increases as 7 increases). 
In fact, the experimenter deals with a probability density function. For example, under 
certain conditions this probability density function may have the form 



x G (xq — A ,997, xq + A ,997), the confidence probability is 0, 997. 
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where a = A 0t997 /3 = A , 6 83 0]- 

Thus, the experimenter knows the state of the system not as a point in a phase space, 
but as a family of confidence intervals with the corresponding confidence probabilities, or, 
equivalently, as a probability distribution in the phase space. One can say that since every 
measurement has an error and the theory of measurement error is based on mathematical 
statistics and probability theory, a description of the state of any classical physical system 
in terms of probability distributions is more natural and direct from the experimental point 
of view. 

Assigning a certain values for the physical quantities, as it is accepted in the traditional 
view on mechanics, is an abstraction and idealization. As every idealization, it works in 
some situations, but in other situations it does not. Of course, we can not neglect the 
measurement errors, when we deal with large times, because the errors grow with the time 
due to the discussed derealization effect. Such large time asymptotics as thermalization 
or recurrence (almost periodicity) are often discussed with the relation to the irreversibility 
problem. So, if we examine the time irreversibility problem, we must consider mechanics in 
its functional, not the traditional, formulation. 

It should be noted that the use of a real-valued probability density function does not 
contradict to the rational-valued measurement results. The probability density function is 
not a directly observable value, rather it is constructed from the directly observable results 
of measurements. See [5] for the detailed description of the construction of the probability 
density function starting from the measurement results. See also [B] for the discussing the 
functional dynamics of a system under often measurements. Since the probability density 
function is not directly observable, it can be real-valued. Like the concept of a material 
point, the concept of a probability density function is also a kind of idealization. But this 
idealization works in a more general case. 

3 Derivation of the Boltzmann equation and the hier- 
archy of space-time structures 

One of the problems of the framework of the irreversibility problem is the derivation of the 
Boltzmann equation from the Liouville equation. Two different derivations have been got 
by Bogolyubov j7] and Lanford [8 J . The main drawbacks of the Bogolyubov's derivation are 
the number of the additional assumptions and the divergencies in the high order corrections 
to the Boltzmann equation [9]. The main drawback of the Lanford's derivation is a small 
time on which the validity of the Boltzmann equation is proved. This time appears to be 
approximately 1/5 of the mean free time [10J. On the contrary, the Boltzmann equation is 
aimed to explain the dynamics of the density function and its convergence to the Maxwell 
distribution on the times much greater than the mean free path. In this section we present 
an alternative derivation of the Boltzmann equation from the Liouville equation. 

Both Bogolyubov and Lanford start with the Cauchy problem for the Liouville equation 




(4) 
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As we have discussed in the previous section, we construct the initial probability density 
function p° from the results of measurements of each particle's position and momentum. 
But practically, when we consider the kinetic events, we do not measure the position and 
momentum of each of iV (which can be of order Avogadro number 10 23 ) particles. In fact, 
there are two space-time scales. The first one is the scale of the particles interaction radius 
r and the particle interaction time r /u (a microscale). Here u is the mean velocity of the 
particles. The second one is the scale of the mean free path I and the mean free time l/u (a 
macroscale or a scale of kinetic phenomena). The Boltzmann equation is valid if r <C I. We 
put p = y — > 0. This parameter in the Hamiltonian (see fl2])) means that particles interact 
on the scale much smaller than the considering scale of I. 

In practice, when we consider the kinetic events, we have a measuring instrument which 
can capture the changes of the functions only on the macroscale. We can construct the 
one-particle probability density function p\{x\) = J u n-i p°dx 2 . . .dxN using this instrument 

(denote also /° = Vp\). But we cannot construct the function p°, because this function 
incorporates the information about the correlations of the particles on the distances of order 
ro which cannot be captured by our instrument. 

So, the Cauchy problem (jl]) has not a direct physical meaning. Instead of the Cauchy 
problem we start with the following problem: 

f = <*■»• 

S? ) M \f2(x 1 ,x 2 ,t- At) - A(zi,t- At)Mx 2 ,t- At)} -> 

f (5) 



(as p — > 0, N — > oo, Np 2 = const, and At — > 0, — = At — > oo) 



/i(x 1 ,0) = / 1 °(xi), 

Hf dx\ . . . dxN < oo, 



k f(x u ...,x N ,t) = f(x l{1) , x i{N) ,t), (i(l), . . . , i(N)) = P(l, . . . , N), 

Here f(x±, . . . , xjv, t) = V N p(xi, . . . , xn, t) (of course, / depends also on p, since p is a 
parameter in the Hamiltonian, the initial function /° also may depend on p); 

f s (xi, t) = V N ~ S / f(x 1} ...,x N ,t) dx s+1 . . . dx N} 



(2) (2) 

s = 1,2; Si is the two-particle Hamiltonian flow, i.e., S t (p{xi,X2) = ip(xi t ,x 2 t), where 
(xit, x 2t) is the phase point (xi, x 2 ) moved along the flow on t] P(l, . . . , N) is a permutation 
of the numbers 1, . . . , N. 

Theorem. Let the function f(xi, . . . ,xn, t) satisfy the problem (TJJ). Let $(r) be mono- 
tonically decreasing function and lim r 7 $(r) = C ^ 0, 7 > 2. Then the function 

r— too 

fi(xi,t) satisfies the Boltzmann equation with the initial function /° in the limit 
p — > 0, — > 00, Np 2 = const = nV (i.e., d ^ converges to the right-hand side of (T3|) 
with fi instead of f in every point x G fly an d t G M). 
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This is the main result of this report. The limit [l — > 0, iV — > oo, iV/i 2 — > const is 
called the Grad or Boltzmann-Grad limit [11] and was used by Lanford. Bogolyubov used 
the thermodinamic limit N — > oo, V — > oo, y = n = const. There is also a derivation of 
Bogolyubov-Boltzmann-type equation for a finite number of particles (for example, for two 
particles), i.e., without the limit iV — > oo, based on the functional approach to mechanics 

Let us discuss the formulation of problem fl5]). The third condition is the initial value for 
the single-particle function fx. The fourth condition means the finiteness of energy, the fifth 
condition means the symmetry of the density function with respect to permutations (i.e., 
the particles are indistinguishable). 

The most interesting and crucial is the second condition. The knowledge of the 
initial single-particle function /i(xi) is not sufficient to get a unique solution for the 
single-particle function fx(xi,t). Since, as we said above, the condition of the form 
f(xi, . . . , xn, 0) = . . . , Xjv) has not a direct physical meaning, we must have some 

additional condition. 

As we said above, there are two scales of space-time in this consideration: the microscale, 
related to the interaction of particles, and the macroscale, related to the kinetic phenomena. 
The two-particle function f 2 relates to the microscopic scale, because it incorporates the 
information about the pairwise correlations of the particles on the distances of order r^. 
The single-particle function fi relates to the macroscopic scale, since it does not incorporate 
the information about particles' interaction and the kinetic theory is expressed in terms of 
this function. The second condition in ([5]) means that the initial value for the microscopic 
function f 2 are assigned by the macroscopic function f\. So, instead of the specification of 
the initial microscopic function f°, we specify only the initial macroscopic function f® and 
impose a condition on the microscopic function: in a certain sense it is subordinated to the 
macroscopic one (in sense that its initial values are assigned from the macroscale). 

Traditionally space-time is defined as a manifold M with a pseudo-Riemannian metrics 
g. In the simplest case this is R 4 with the usual Lorentz metrics rj (the Minkowski space). 
However, when we speak about a space-time, we keep in mind some characteristic length and 
time. For example, in continuum mechanics the density is defined as d ™y^ = lim , 

where Am(r) is the mass of the volume AV near the point r G R 3 . In spite of the mathe- 
matical limit AV 0, i.e., the limit of arbitrarily small volumes, we keep in mind that AV^ 
is a macroscopic volume, i.e., it contains a large number of atoms or molecules (otherwise 
the definition is not valid). 

So, always when we deal with some space-time, we keep in mind not only the space-time 
(R 4 , n) itself, but also a characteristic length L, which tell us about the averaging scale 
(the characteristic time is given by L/u, where u is the mean velocity). In our example we 
average over the molecular structure of the medium. Exactly this allows us to perform the 
limit A V — > 0. The triple (R 4 , n, L) we will call the space-time structure. 

In our problem of the derivation of the Boltzmann equation we also actually deal with 
two distinct space-time structures. The first one is the space-time structure (G x R, n, I) of 
macroscopic (kinetic) phenomena with the coordinates (q, t) and the characteristic length I 
(the mean free path). The second one is the space-time (— x R, n,r ) of microscopic phe- 
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nomena, like particle interactions, with the coordinates (£, r) and the characteristic length 
I (the particle interaction radius). These space-time structures are related to each other by 
the scale transformation 

£=-, T=-, fi=-r->0. 6 
/i /i I 

With such a transformation infinitesimal region of the macroscopic space-time may corre- 
spond to the infinitely large region of the microscopic space-time. Exactly this situation 
we can see in the second condition in ([51): At — > 0, but At — — — >• oo. Of course, these 
two scales meet each other in the collision integral of the Boltzmann equation: a collision 
is considered as a point and momentary act on the macroscale, but it takes place on the 
infinite space during the infinite time on the microscale. 

The great disparity between microscopic and macroscopic scales as one of the origins of 
the irreversible macroscopic behaviour was pointed out in [121 [TJ]. The used rescaling of 
the space-time ([6]) is typical for the derivation of the Boltzmann equation from the Liouville 
equation [TC)| [15] . In the case of lattice dynamics this rescaling was used in [161 E] for the 
derivation of kinetic and hydrodynamic-type equations. Our proposition is to introduce the 
subordination of different space-time structures expressed in the form of the second condition 
in (0), which gives a new way of derivation of the Boltzmann equation. 

One can say about a kind of hierarchy of the space-time structures: the initial values 
for the processes on the microscopic space-time structure (interactions of the particles) are 
assigned from the processes on the macroscopic space-time structure (kinetic phenomena). 
Note that the idea of the hierarchy of times (namely, the microscopic, kinetic and hydrody- 
namic times) in a slightly different sense was first proposed by Bogolyubov [7]. 

Remark. We can note that in our limiting process the overall mass of the gas M = mN tends 
to infinity, since m is constant and N tends to infinity. However, we can also rescale the 



i 2 r 

way $ M (r) = yU 2 $(r) as it was proposed by Grad [llj (we substitute m by m M and $( 

by $ M ( ^'~ g ^ ) in fl2])). The Hamiltonian system with these additional rescalings can be easily 
reduced to the considered one. So, the gas is considered as an infinite number of negligible 
mass particles in the limit. This exactly corresponds to the macroscopic intuition. 

Finally, we would like to note that the functional approach can be useful in the mothod of 
molecular dynamics simulations. At present the molecular dynamics follows the Newtonian 
approach and simulates the movement of material points. On this way it is hard to obtain 
the properties of complex and, moreover, biological systems, which constitute the aim of 
the molecular dynamics. Furthermore, the problem of uncontrolled cumulative errors in 
numerical integration is known. The simulation of the motion in terms of the probability 
density function seems to be more appropriate for obtaining the properties of the complex 
and biological systems and more justified from the computational point of view. Also the 
present report suggests that the initial conditions for the probability density function should 
be chosen not in an arbitrary way, but rather in the form suggested by the system as a whole. 
In our case such conditions are given by (E]). Let us note that these conditions are statistical 
in essence, they are formulated for the probability distribution and cannot be reduced to the 
conditions on the initial positions and momenta of the individual particles. 



mass of a single particle in the way m M = ji m and rescale the interaction potential in the 



8 



4 Conclusions 



In this report we have reviewed the irreversibility problem and the functional formulation of 
mechanics. This formulation is based on the concept of a probability density function. The 
fundamental equation of this formulation is the Liouville equation for the probability density 
function. The functional approach to mechanics is justified by the measurement theory. We 
have proposed a derivation of the Boltzmann equation from the Liouville equation based on 
some ideas of the functional formulation and the measurement theory. The main features 
of this derivation is distinguishing of two scales of space-time (micro- and macroscopic or 
kinetic) and subordination of the processes on the microscopic space-time structure to the 
processes on the macroscopic one. 

According to the traditional paradigm of the mathematical physics, the dynamics is 
completely determined if we know the law of motion, i.e., a differential equation, and the 
initial values for it. However, the initial values themselfs are understood as something 
external to the equations of mathematical physics ("As regards the present state of the 
world. . .the laws of nature are entirely silent" [H]). We propose another picture: the initial 
values for a given level of nature are assigned from the higher level. 

So, instead of the reductionism, which claims the reducibility of all levels of the nature to 
the most microscopic level, we propose a hierarchical picture of the world: the lower levels 
of the nature are subordinated to the higher ones. 
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